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ABSTRACT 


Silicon-on-insulator bipolar transistors fabricated using the Harris UHF-1 
process, were irradiated at room temperature with 30 and 60 MeV electron 
beams. Some of the transistors on each die were configured and biased as a 
simple operational amplifier (opamp), one was placed in a common emitter type 
circuit and the remaining were biased to measure transistor parameter 
degradation. The purpose of this setup was to observe the total dose effects of 
the transistor and of an opamp on the same die in order to derive a more accurate 
model of an opamp under total dose conditions. This investigation was 
successful in conducting in-situ measurements of opamp gain and 3dB jfrequency 
while also measuring the current gain of similar transistors on the same die. 
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I. INTRODUCTION 


The subject of radiation effects on electronic material 
has been a heavily studied subject since the early 1960's and 
especially since the Starfish nuclear event (1962) which 
increased the radiation levels in the Van Allen belts, 
resulting in the failure of several satellites in orbit [Ref. 
1, p. 1527]. The degradation of bipolar transistor current 
gain caused by ionizing radiation was the cause of these 
satellite failures [Ref. 2, p.l540] . Interest in this subject 
is still high and rightly so, since almost every aspect of our 
lives depends on electronic devices, more of which are being 
deployed on spacecraft. 

Now there is even a greater variety of electronic 
components and circuits with an ever increasing complexity and 
range of failure mechanisms. These circuits and devices are 
being used extensively in space, where there is often a harsh 
radiation environment, even in the absence of a nuclear 
detonation. In order to ensure that spacecraft can complete 
their mission safely and effectively, there must be a thorough 
understanding of the mechanisms that cause failure in 
electronic components. 

Regardless what is sent into space, it will be exposed to 
a certain amount of radiation over its useful life. A 
designer of space products needs to know how materials, 
specifically semiconductors and shielding materials, will 
react given a certain dose of radiation. In order to 
fabricate an effective device to withstand radiation, 
radiation models need to be developed for the particular kinds 
of devices utilized. 

In order to design a device or circuit for in space or 
other radiation environment there must be an awareness of the 
types and levels of radiation in the operating environment. 
The nature of the space radiation environment (the environment 
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of concern in this paper) cannot be understood with a simple 
model. For example: the total dose received in ten hours 
inside Spacelab on a low altitude Space Shuttle flight was 10'^ 
rad while it was expected to receive up to 10® rads in one 
hour at the closest approach to Jupiter on the Galileo flight 
[Ref. 1, p.l530] . Orbits in the somewhat homogeneous trapped 
radiation belts can experience sudden changes in dose due to 
solar flares occurring during various times of the year. 
Figure (1.1) shows the dose/year for aluminum as a function of 
altitude for a particular orbit. As can be seen, the nature 
of space radiation will lead to the discovery that it is a 
function of both space and time. 



Figure 1.1. Dose rate versus altitude for a circular orbit 
at 600 inclination. From Ref. [2, p. 1511] . 

The space radiation environment for orbits up to the 
geosynchronous orbits can be broken down into three sources: 
trapped radiation belts (Van Allen belts) consisting of 
protons and electrons, solar flares and cosmic rays (mostly 
from outer galactic sources). These sources can result in 
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dielectric ionization (total dose effects), transient 
ionization and displacement damage. The nature of the space 
radiation environment is further discussed in Chapter II. 

Of the three above effects, displacement damage is most 
predictable. This affects primarily solar cells and wide-base 
bipolar transistors (used mostly in power applications). 
Transient ionization (single particle effects) causes what is 
known as single event upset (SEU) in MOS type memory devices. 

Dielectric ionization or total dose effects were the 
first to cause failure in the space radiation environment and 
continue to cause problems in more advanced and high density 
microelectronic circuits. In this research, where bipolar 
devices are of concern, the major effect of ionization is an 
increase in surface recombination velocity which ultimately 
degrades current gain. Chapter III discusses this effect in 
further detail. 

From the above discussion, it can be seen that radiation 
testing is a necessity for the design of space applications. 
The goals of radiation testing are best summarized by Pease: 

The objectives of radiation testing are 
threefold: 1) understand the mechanisms of the 
interaction of the radiation with electronic 
materials and how these effects relate to device 
failure; 2) characterize the response of specific 
device types and technologies for use in part 
selection for specific systems application; and 3) 
determine the acceptability of production lots for 
use in deployable systems. [Ref.2, p. 1514] 

The goal of this research is reflected by the first objective 
above. More specifically, we want to form an accurate model 
for a given technology that can predict the response of 
operational amplifiers to total dose radiation. The approach 
taken is to first investigate the response of individual 
components of the opamp (transistors in this research) and 
apply these results to form a model for an opamp by way of 
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superposition of the responses of individual components. It 
is imperitive then, to validate this data by exposing the 
opamp to the same radiation environment as the individual 
components. To accurately do this, the opamp and separate 
components must be on the same die and certain device 
parmeters and circuit design of the opamp must be known (this 
can sometimes be difficult due to manufacturer proprietary 
reasons). This research is primary directed to a proof-of- 
concept of the above. 

The obvious method for radiation testing would be to put 
the device or circuit in the actual environment and monitor 
it. This is done in fact (Combined Release and Radiation 
Effects Satellite (CRRES)), however the time required and 
costs can be excessive and the type of data that can be 
acquired is limited [Ref. 2, p.l514]. Therefore, alternative 
solutions must be sought. 

Reproduction of the space environment in the laboratory 
can be difficult and expensive. More commonly, a single 
radiation source is used to isolate the dominant effect in a 
particular space environment. The damage in a total dose 
environment is mainly caused by ionization energy absorbed by 
the material. This type of damage can be simulated with 
electrons or photons. More costly proton sources can also be 
used, but the disadvantage is that displacement effects must 
be differentiated from the ionization effects. Table 1.1 
summarizes different total dose simulator characteristics and 
their comparison with the space environment. 

Since this research is concentrating on total dose 
effects, the simulation used in this experiment was conducted 
using 3 0 and 60 MeV electrons generated by the NFS Linear 
Accelerator (LINAC). In addition, in-situ (sometimes referred 
to as in-flux) testing was utilized. This method involves 
taking the device parameter measurements while it is being 
irradiated under normal operating conditions. Another method. 
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Energy 

rad (Si 02 )/s 

Width 

Dynamitron 

electrons 

2.5 MeV 
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5-60 MeV 
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10 keV to 
> 100 keV 

10’-10^ 

dc 

Space 

electrons. 

< 1 MeV to 

<10"’ ave. 

variable 

photons, 

> 100 MeV 




protons 





Table 1.1. Total dose simulator 
characteristics compared to the space 
ionization environment. From Ref. [2, p. 
1516] 


step-stress testing, involves taking measurements, irradiating 
the device and then retaking the measurements for comparison. 
This method is a lot simpler than the in-situ method but the 
results are considered to be less precise [Ref. 2, p. 1515]. 

The testing in this research consisted of irradiating 
silicon-on-insulator (SOI) bipolar transistors biased for 
current gain measurements, implemented in a simple operational 
amplifier configuration and in a common emitter configuration. 
All transistors were co-located on a single die to ensure the 
same radiation exposure. The complete description of this 
experiment is contained in Chapter IV. In Chapter V the 
results are discussed and the conclusions and recommendations 
for further research are found in Chapter VI. 
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II. RADIATION ENVIRONMENT AND ITS EFFECTS 


A. FORMS OF RADIATION AND THEIR INTERACTION 

In the space radiation environment there are many forms 
of radiation that may affect vital electronic components in 
spacecraft (specifically satellites). Electrons, protons, 
gamma rays, alpha and beta particles and heavy ions are of 
prima 3 ry concern. These forms of radiation can be categorized 
into two types: photons (electromagnetic waves) and 
particles. 

Electromagnetic waves consist of discrete quanta of 
energy hv, where h is Planck's constant and v is the 
frequency. These quanta are called photons and are used to 
describe the particle-like nature of electromagnetic 
radiation. A photon can be considered a particle with zero 
rest mass and neutral charge. Two examples of photons are 
gamma rays and X-rays. There are three types of interactions 
a photon can have on matter: the photoelectric effect, 
Compton scattering and pair production. 

The photoelectric effect was first observed by Heinrich 
Hertz in 1887 and further investigated by Lenard in 1900. 
When light of frequency v is incident on a clean metal 
surface, electrons are emitted with a maximum energy described 
by the Einstein equation: 

E = e Vq = = hv - (2.1) 

where V© is the stopping potential, m is the mass of an 
electron, v is the velocity of the electron and 4)o is the 

energy required to remove an electron from the surface of the 
metal (work function). Figure (2.1) shows the experimental 
setup for the photoelectric effect and a plot of the observed 
photocurrent vs. stopping potential Vq . [Ref. 3, pp. 107-109] 
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Light 



Figure 2.1. Photoelectric effect. From Ref. [4, p. 66]. 

From Equation (2.1), it can be seen that for the 
photoelectric effect to occur, the incident photon must have 
energy at least greater than some threshold energy hVo, 

regardless of the incident photon intensity. In this case, 
the number of emitted electrons (sometimes called 
photoelectrons) are proportional to the incident light 
intensity and have frequency v-Vq. Otherwise the incident 

light energy will just increase the kinetic energy of the 
electrons, resulting in a temperature increase without the 
release of electrons. Note that Equation (2.1) only gives the 
maximum energy given the incident monochromatic frequency. 
The actual energy of the electrons leaving the surface may be 
less due to energy loss transversing the finite thickness of 
the metal. [Ref. 3, pp. 107-108] [Ref. 4, p. 66] 

Compton scattering is another phenomenon demonstrating 
the particle nature of electromagnetic radiation. When a 
photon of energy h\ and momentum p=hv (X=c/v) collides with 
an electron, the photon is scattered by an angle of 0 and 
loses part of its energy to the recoiling electron. This 
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transfer of energy is reflected by a decrease in the frequency 
of the photon and is related by the following equation: 

= (1-COS0) (2.2) 

where and X^ are the photon wavelengths before and after 
the collision respectively. Figure (2.2) illustrates this 
effect. Ea=mc^ is the energy of the electron at rest and 

E={eI+pIc^)^^^ is the energy of the electron after the 
collision. [Ref. 4, pp.115-116] 



Figure 2.2. Compton scattering. From Ref. [3, p. 116]. 

When a photon of sufficiently high energy passes nearby 
an atomic nucleus, the nucleus recoils and an electron- 
positron (has properties of an electron but possesses a 
positive charge) pair is formed. This process is called pair 
production or materialization. The positron is only 
transitory and will soon interact with an electron forming two 
or three gamma rays (the formation of one gamma ray would 
prohibit the simultaneous conservation of energy and momentum) 
[Ref. 5, p. 922]. The incident photon is annihilated in the 
process. If the recoil of the nucleus is absorbed by an 
electron, a triplet is formed consisting of two electrons and 
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a positron. Figure (2.3) is a photograph of hydrogen bubble 
chamber used for observing elementary particles. In upper 
half of Figure (2.3), three tracks can be seen to originate 
from a common point. This is where the photon interacts with 
the nucleus and is annihilated producing a triplet consisting 
of an electron-positron pair (the circular paths) and another 
electron. In the lower half of the figure, two tracks can be 
seen, originating from another common point (the circular 
paths are not visible). This is the electron-positron pair 
without the extra electron. Observe that the triplet 
electron-positron pair tracks curve in opposite directions. 
Since the bubble chamber is in a magnetic field, comparison of 
the two tracks indicate that the two particles have opposite 
charges. Both particles are shown to ionize like electrons in 
the chamber, thus proving the existence of a positron. The 
threshold energy for this process is 1.022 MeV. [Ref. 5, p. 
589] [Ref. 6, p. 55] [Ref. 7, p.479] 

The energy of the photon determines which of the above 
interactions dominate. In the case of silicon, the 
photoelectric effect dominates at photon energies less than 50 
keV. Between 50keV and 20MeV Compton scattering dominates. 
When photon energies are greater than 20MeV, pair production 
dominates. [Ref. 8, p. 1444] 

Charged particles interact with matter quite 
differently than photons. Elastic collisions occur between 
charged particles and atoms called Rutherford scattering. In 
particular, the nucleus of one atom will scatter particles at 
a certain angle 0 that are incident at a distance b from the 
axis of the atom (see Figure 2.2). This occurs only if the 
particles do not have enough energy to penetrate the nucleus 
of the atom. 
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Figure 2.3. Pair production in a hydrogen 
bubble chamber. Tracks show formation of 
triplet and pair. From Ref [6, p. 56]. 
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The following Rutherford scattering equation describes this 
interaction [Ref:2, p.23]: 

)2 sinBdB 

A%ejnv^ sin^( —) (2.3) 

Note that Equation (2.3) is derived assuming that the nucleus 
is a point charge. This is not the physical case but the 
equation holds true so long as the incoming particle does not 
penetrate the nucleus. [Ref. 3, p.138-139] 



Figure 2.4. Rutherford scattering. From Ref. [1, p 

139] . ' F 

B. RADIATION IN SPACE 

1. Historical Background 

In the 1950's, the scientific community's perception of 
the earth's magnetic field was that it was like a bar magnet 
as in Figure (2.5) and that it became weaker as it extended 
thousands of miles into space. 
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In addition, cosmic rays were known to reach the earth at 
the time but were thought to originate from somewhere in space 
other than the sun (which was believed only to emit visible 
light) since the rays were isotropic. It wasn't until the use 
of rockets and orbiting satellites that the nature of space 
radiation was revealed. [Ref. 9, p.2] 



Figure 2.5. Early view of the 
Earth's magnetic field Ref 
[13, p. 4-2] . 


After much political tuirmoil with the U.S. satellite 
program Project Vanguard, Explorer I was finally launched on 
January 31, 1958 carrying a Geiger counter provided by 
Professor Van Allen's group in Iowa. There was no recorder on 
board Explorer I, so data had to be recorded by receiving 
stations during satellite overpasses which were approximately 
2 minutes in length. The data received from Explorer I 
revealed that over South America at high altitudes of over 
2000 km, there appeared to be no cosmic rays. However, at low 
altitudes of 500 km, the normal cosmic ray count of 30 
counts/sec was observed. The Iowa group could not come up 
with an explanation for this. It wasn't until Explorer III 
was launched on March 26, 1958 that they found an answer. 
Since this satellite had a recorder on board, continuous data 
could be taken. At the beginning of the pass over South 
America, the count rate rapidly increased from the normal 
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cosmic ray level to recorder's limit of 128 counts/sec and 
then dropped to almost zero after seven minutes. About 10 
minutes later, it increased sharply to 128 counts/sec. Carl 
Mcllwain of the Iowa group explained that the counter had 
reached an area of a very high particle density and that the 
"dead-time effects of the counter reduced the count rate 
essentially to zero" [Ref. 9, pp. 5-9]. This explained the 
apparent absence of cosmic rays. In addition. Van Allen 
theorized that the increase in radiation over the cosmic level 
was due to trapped charged particles. This marked the 
beginning of the formulation of the Van Allen radiation belt 
model. Figure (2.6) shows Van Allen's first map of the 
radiation belt derived from two later satellites. Explorer IV 
and Pioneer III, which was published in 1959 [Ref 10, p. 39] . 



Earth 

radii 


Figure 2.6. Van Allen's original map of the 
radiation belt. Constructed with data from 
Explorer IV and Pioneer III. From Ref [8, p. 39] . 


14 



2. Types of Space Radiation and Its Sources 

Radiation encountered in space can be broken down into 
three types or "sources": cosmic rays, geomagnetically trapped 
radiation near the earth (Van Allen belts) and solar plasma 
[Ref. 10, p.l] . Spacecraft will encounter some or all of 
these types of radiation depending on its orbit. 
Understanding of the spatial distribution and effects of each 
type of radiation is essential when designing systems to be 
utilized by the spacecraft. 

a. Cosmic Rays 

Cosmic rays were first observed by Austrian 
physicist Victor Hess in 1911, using manned balloons at a 
maximum altitude of 17,500 feet. Since then many others 
including Robert Millikan conducted similar measurements at 
even higher altitudes, confirming "an extra terrestrial source 
of penetrating radiation." [Ref. 11, p.7] 

A small amount of cosmic rays originate from the 
sun, where there are minor changes in flux from day to night, 
but the majority are from the Milky Way galaxy and even beyond 
where the sources are not exactly identified. Cosmic rays 
consist mostly (about 90%) of high-speed protons (the nuclei 
of the hydrogen atom), helium nuclei or alpha particles (9%), 
electrons (less than 2%) and gamma rays. [Ref. 7, p. 461] 

Cosmic ray protons incident on the earth's 
atmosphere collide with nitrogen and oxygen nuclei forming 
secondary particles called TT-mesons (mesons) which then decay 
rapidly (about 10‘®sec) to ^-mesons (muons). The ^-mesons 
decay in about 2xl0’®sec to electrons (energies up to 50 MeV) 
and uncharged ^-mesons decay into two gamma rays within 10 
sec. Hence, inbound cosmic rays to the earth's atmosphere 
form an avalanche or showers of electrons, gamma rays and 
various short-lived elementary particles. Figure (2.7) 
ilustrates this effect. [Ref. 9, p. 70] [Ref. 11 pp. 12-13] 
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Figure 2.7. Cosmic Ray- 

Avalanche. From Ref. [10, p. 

5] . 

Neutrons are also produced by these cosmic ray proton 
collisions with the atmospheric nitrogen and oxygen nucleus. 
The neutron is radioactive with a mean life of about 1000 sec 
and decays as follows: 

n - p + p- + Vg (2.4) 

The antineutrino was first introduce by Pauli 

in 1930 to explain the absence of conservation of energy in p 
decay and was observed in the laboratory in 1957. As far as 
space radiation is concerned, the important part of this decay 
is the resulting proton and p decay which is essentially the 
emission of a radioactive electron with a half life of about 
10.8 min. The energy spectrum of the emitted electron is 
independent of the neutron's kinetic energy (assuming neutron 
energies less than lOMeV) and will have a p decay spectrum as 
shown in Figure (2.8). The resulting proton, on the other 
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hand will have kinetic energy almost identical to the parent 
neutron. This is because the electron's kinetic energy is 
almost the same as the mass difference between the neutron and 
electron. [Ref. 9, pp. 71-74] 



Electron energy, kev 

Figure 2.8. Beta decay spectrum. 

From Ref. [9, p.74]. 

Neutrons are also produced by the sun predominately 
by knock-on reactions in hydrogen and helium. However, it 
should be noted from the above discussion that very few 
neutrons actually make it to the earth's surface due to the 
short lifetime of the neutron. Figure (2.9) shows the 
measured neutron flux at various altitudes for various 
energies. [Ref. 9, pp. 72-73] 

Jb. Solar Plasma 

Plasmas are regions of particles, charged and 
uncharged, that when exposed to a static electric fields 
redistribute themselves so that they are shielded from the 
field. On the boundary of plasmas are strong field regions 
called sheaths. An important property of plasmas is that they 
tend to remain electrically neutral. [Ref. 12, p.291] 

Solar plasma, as the name implies, originates from 
the sun and consists of equal numbers of electrons and 
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Figure 2.9, Neutron-energy spectrum at different distances 
above the geomagnetic equator. The top of the atmosphere is 
labeled 0 Rg (approx. 100km) . From Ref. [7] . 


positive ions. It is actually a continuation of coronal 
expansion and begins at a very slow speed at the corona but 
becomes supersonic at 10-20 solar radii (6.96x10^ km) . This 
is what comprises solar wind whose effects will be described 
in the next section. 

c. Van Allen Belts 

Strictly speaking the radiation belts are not 
actually a source of radiation. The effect and interaction 
that the magnetosphere and atmosphere has on incident 
radiation can give the appearance that the radiation belts are 
a source. Understanding these radiation belts is critical 
because of the number of satellites in this region. Low earth 
orbits (LEO'S) include altitudes of 90 to 500 miles. 
Satellites in LEO'S include weather satellites, space shuttle, 
space telescope, earth resources satellites and space 
stations. In addition. Global Positioning Satellites (GPS) 
orbit at about 12,000 miles, which is within the outer Van 
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Allen belt. Because of its complexity and importance. Van 
Allen belts are treated separately in the next section. 

3. The Magnetosphere and the Van Allen Belts 

The natural trapped radiation environment in orbits about 
the earth (altitudes to about 32,000 km) as we know it today 
consist of electrons protons and a lesser number of heavy, low 
energy ions. The earth's magnetic lines of force are 
constrained to a region of influence or cavity called the 
magnetosphere by the continuous flow of charged particles (the 
solar plasma) from the sun called solar wind. The outer 
boundary of this region is known as the magnetopause (Refer to 
Figure (2.10)). The pressure of solar wind shapes this 
magnetosphere. As a result, the daylit side is compressed and 
the night side is expanded as shown in Figure (2.10). 



Figure 2.10. Magnetosphere. From Ref. [14, p. 1424 ]. 
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The magnetosphere can be divided into five regions 
according to the particles predominating as shown in Figure 
(2.11) . The boundaries of these regions are not to be 
interpreted as sharp interfaces but diffused areas changing in 
position due to various solar and other variations. 


REGION 

SOLAR FLARE 
PROTONS 

TRAPPED 
PROTONS 
OUTER ZONE 
ELECTRONS 

INNER ZONE 
ELECTRONS 

1 2 3 4 S 6 7 8 9 10 11 12 13 14 

L (EARTH RADII) 

R = L cos^A 
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TRAPPED 
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SOLAR FLARE 
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Figure 2.11. Distribution of charged particles in the 
magnetosphere. From Ref. [8, p. 1425]. 

The path of charged particles from the sun or other 
sources approaching the magnetospere will be altered by the 
field and tend to follow the field's lines of force. 
Particles below a critical energy become trapped and tend to 
spiral around these lines and bounce back and forth between 
reflection or mirror points (regions of maximum magnetic field 
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strength for a given path). Figure (2.12) illustrates this 
motion of particles. The charged particles (electrons and 
protons) that have become trapped within the plasmasphere 
region form what are called the Van Allen belts. 



Trajectory of 
trapped particle 




Mirror point ^ 

(Pitch angle of helical trajectory=90 ) 


Magnetic field line 





Figure 2.12. The motion of particles inside the trapped 
radiation belts. From Ref. [8, p. 1425]. 

A useful coordinate system, introduced by Mcllain in 
1961, relates the magnetic field B to an adiabatic invariant 
(constant only for very slow changes of the variables 
involved) called the integral invariant J. The integral 
invariant is obtained from the longitudinal invariant J [Ref. 
9, p.49] . The quantity J is the length of the field line 

between the reflection points described above. This 
coordinates system was designed to relate measurements made at 
physically different locations. Mcllwain introduced theL 
parameter which is a dimensionless ration of the earth's 
surface, approximately the geocentric distance of a field line 
in the geomagnetic equator. The L parameter is given by: 

L = i?(cosA) .5) 

where A is the magnetic latitude [Ref. 13, p. 5-8]. 
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The transformation of B is now given by: 


B = 


M 


3R 

L 


( 2 . 6 ) 


where M is the geomagnetic dipole moment [Ref. 13, p. 5-8]. 

There is a particle flux dependence on latitude and 
altitude. In the lower latitudes, incoming particles are 
deflected sharply to follow the magnetic lines of force. 
Hence, only very high energy charged particles reach the low 
altitudes at low latitudes. At high latitudes, particles are 
already travelling along the lines of force (or nearly so) and 
therefore little deflection takes place. Even low energy 
particles achieve low altitudes at high latitudes. This is the 
reason for the depression of the Van Allen belts toward the 
earth's magnetic poles. Generally, the least radiative 
environment is at low inclinations (below 28°) and low 
altitudes. This would correspond with the Low Earth Orbits 
(LEO). Figure (2.13) summarizes the radiation doses at 28°. 
[Ref. 14, p.1424] 

4. Outside the Van Allen Belts 

Geosynchronous equatorial orbits (GEO) at an altitude of 
22,300 nautical miles, are essentially the highest orbits used 
by satellites. Many critical satellites such as military and 
commercial communication satellites. Defense Support Program 
(DSP) and Tracking and Data Relay Satellite System (TDRSS), 
utilize these orbits. Thus, it is essential to understand the 
radiation environment in this region of space. [Ref. 16, p. 
141] 

The outer limit of the trapped radiation would have to be 
where particles would no longer bounce between mirror points 
as was shown in Figure (2.9). This distance is about 10-12 
earth radii on quiet days and does not change much with 
changes in the solar wind. However there are some regions in 
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Figure 2.13. Summary of Radiation Doses. From Ref. [15, p. 
2.2-4] . 

which particles whose path starts in the day side and overlaps 
in the night region where the field is no longer closed. 
These particles will be lost in the geomagnetic tail. This is 
called psuedotrapping and the regions where this occurs in 
shown in Figure (2.14). [Ref. 11, pp. 322-326] 

The plasma sheet is made up of hot plasma that contains 
equal parts of hydrogen and oxygen ions. The electrons in 
this sheet have energies on the order of lOkeV and can reach 
as high as BO-lOOkeV during magnetic substorms. The plasma 
sheet is the normal environment for geosynchronous orbits. 
[Ref. 13, p. 5-17] 
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Figure 2.14. Geomagnetic-tail in the noon-midnight plane. 
From Ref [9, p. 334] . 
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III. TRANSIENT EFFECTS OF RADIATION ON ELECTRONIC SYSTEMS 
(TREE) 

A. INTRODUCTION 

Transient Effects of Radiation on Electronic Systems 
(TREE) is a general term associated with the effects of total 
dose, dose rate and neutron displacement in electronic 
devices. This chapter will attempt to provide a very brief 
introduction to the concepts of semiconductor physics and 
apply them to the basic operation of a bipolar transistor. 
Then the Transient Effects of Radiation on Electronic Systems 
(TREE) will be discussed for the specific case of the bipolar 
transistor. 

B. OVERVIEW OF SEMICONDUCTOR PHYSICS 

Silicon, a group IV element in the periodic table has 
four valence electrons which form covalent bonds with 
neighboring atoms resulting in a diamond crystalline structure 
called interlocking-face centered cubic [Ref. 17, pp. 19-20] . 
Figure (3.1) illustrates the geometry. Figure (3.2) gives a 
2-dimensional view of the atom array showing the sharing of 
electrons in the covalent bonding. This is called intrinsic 
silicon (no impurities) and the nature of the bonding results 
in approximately 5x10^^ atoms/cm^ [Ref. 19, p. 36-1]. 

The valence electrons can be in either of two energy 
states, the conduction band or valence band. An electron in 
the valence band that is excited by some outside energy source 
can jump to the conduction band, leaving a vacancy called a 
hole which results in a positive charge in the bond. Another 
electron in the valence band can fill this hole, creating 
another hole in its original position. This repeated process 
can generate a net movement of electrons in one direction and 
holes in the opposite direction. This process is called 
electron-hole generation. In an intrinsic crystal at thermal 
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Figure 3.1. Diamond crystal structure. 

Shaded atoms are of the original face- 
centered cube. From Ref. [18, p. 8]. 

equilibrium, the number of holes in the valance band equal the 
number of electrons in the conduction band. [Ref. 20, p. 9] 
Since electrons and holes obey Fermi-Dirac statistics 
according to quantum theory, the probability that an electron 
will be in a certain energy state is given by the Fermi-Dirac 
distribution function: 

f(F) = - _ 

where f{E) is the probability the electron will be 
state E, Ef is the Fermi energy level at which the 

has a probability of M, k is the Boltzman constant 
eV/°K) and T is temperature in °K. The density of 
states per unit volume or number density of states 


(3.1) 

in energy 
electron 

(8.63x10-^ 
allowable 
is: 
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Figure 3.2. Two dimensional version of 
the silicon lattice illustrating the 
covalent bonds. The shaded circles 
represent the Si nucleus the inner two 
orbitals. The solid circles represent the 
valence orbitals. Ref. [18, p. 16]. 


N,{E) = - 


(3.2) 


where m* is the density of states effective mass of an 


electron in the conduction band and E^ is the energy at the 

bottom of the conduction band. The effective mass is an 
adjustment to the carrier's mass to take into account that the 
carrier not only experiences force from an external electric 
field, but also forces from other electrons and atoms within 
the material. The effective mass is not only material 
dependent, but also energy band dependent. [Ref. 19, p. 36-3] 
[Ref. 21, pp. 10-11] 













The actual number density of states is given by: 

n^= r N^{E)f{E)dE (3.3) 

The upper limit is taken to infinity vice the energy of the 
upper conduction band because the Fermi-Dirac distribution 
function contributes little in the upper limit. Normally 
E^-Ef>3kT so that the integral can be evaluated to: 

= 2 {2%mlkT/h^) exp{Ef-E^)/kT (3.4) 

The result is identical for holes, except that the effective 
mass for the electrons in the conduction band is replaced for 
the effective mass of the hole in the valence band and E^ is 

replaced with E^, the energy at the top of the valance band: 

= 2 {2nmpkT/h^) exp{Ef-E^) /kT (3.5) 

The only carriers in an intrinsic semiconductor are electrons 
which are equal to the number of holes. Therefore no=Po in an 
intrinsic semiconductor and their product is a constant for a 
given material at a given temperature. The intrinsic carrier 
density is given by: 

Po 

The intrinsic carrier density of silicon at room temperature 
is approximately 1.4X10^° carriers/cm^. [Ref. 21, pp. 11-12] 
The carrier density in intrinsic semiconductors is not 
high enough for the currents required for proper operation. 
To increase the carrier density, intrinsic semiconductors are 
doped (diffused) with impurities. Usually a group V element 
such as phosphorus or arsenic is used to produce n-type 
material. Group V element atoms are called donor atoms 
because four of the five outer electrons form the covalent 
bonds with the silicon and the fifth electron can be easily 
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ionized from the donor atom which will add to the intrinsic 
carrier density. The donor (and acceptor) density, Nj, can 
vary from 10^^ atoms/cm^ to a maximum of about 10^° atoms/cm^ 
for silicon. This means that for in the lattice structure of 
the silicon crystal, one silicon atom out of every 500 to 10^° 
is replaced by the donor atom. [Ref. 21, p.l2] [Ref. 20, p. 
10 ] 

For p-type material, silicon is doped with group III 
acceptor atoms such as aluminum or boron. The outer shell of 
these atoms contain three electrons, therefore a hole will be 
produced when these acceptor atoms removes an electron from 
the valence band to form the covalent bond in the crystal 
structure. [Ref. 20, p. 10] 

New energy levels are formed just inside the forbidden 
gap by the introduction of impurities in the semiconductor. 
In n-type material this energy level is just below the 
conduction band and for p-type it is just above the valence 
band. These new energy levels are shown in Figure (3.3) . 



Intrinsic rr-T ype P~Type 


Figure 3.3. Illustration of the shifting 
of the Fermi energy level due to dopants. 
From Ref. [21, p. 14] . 


The existence of these new energy levels bring about a change 
in the Fermi energy level which is given by: 
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Wb - ^a\ 


(3.7) 


Ef = Ef. ± kT In 

where Ef^ is the intrinsic fermi level given by Equations 
(3.4) and (3.5) and and are the donor and acceptor 

concentrations respectively. It should be noted that since 
the doped semiconductor is neutral electrically, the following 
is true: 


n + Nj^ = p + (3.8) 

where n and p are the electron and hole carriers 
respectively. Equation (3.7) makes use of this equation and 
the assumption that for n-type material. [Ref. 21, 

pp. 14-15] 

At 3 0 0°K the majority of the impurity atoms are ionized, 
meaning that in n-type material the impurity atoms have given 
up their electrons and in the p-type material the impurity 
atoms have taken the electrons from the valence band to 
complete their covalent bonds in the crystal structure. Since 
most of the carriers in the doped material are generated by 
the impurity atoms, the majority carriers (at room 
temperature) in n-type material are electrons and in p-type 
material the majority carriers are holes. It should be noted 
that the impurity ions in both types of semiconductor material 
considered to be immobile. Additionally it is assumed 
that most of the carriers are due to the donor and acceptor 
atoms and very little contribution is from the intrinsic 
semiconductor atoms. Table 3.1 summarizes the makeup of the 
charged particles in an extrinsic semiconductor. [Ref. 22, p. 
325] 

From the above discussion, it should be clear thatnsi\7p 
and p^N^ for n-type and p-type material, respectively. Using 
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/f-type 

p-type 

Impurity (dopant) 

Impurity concentration (cm''-’) 
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(No) 
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(Na) 
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Electron 

(n«) 
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(Pp) 

Minority carrier 

Minority carrier concentration (cm“’) 

Hole 

(p») 

Electron 

(ftp) 

Minority carrier concentration 
at thermal equilibrium (cm”’) 

(Pna) 

("pa) 


Table 3.1. Summary of carriers and 

impurities. From Ref. [20, p. 11]. 


this and Equations (3.4) through (3.8), the following can be 
derived for n-type and p-type material: 


n s 




(3.9) 


From both Equations (3.8) and (3.9), it can be seen that 


iii = np 


(3.10) 


also holds for extrinsic materials. Equation (3.9) is not 
valid for extremely high doping concentrations that are over 
10^®/cm’^ or for extreme temperatures. [Ref. 22, p. 335}] 

When an electric field E is applied to a semiconductor, 
the average velocity of an electron, v^, and the average 
velocity of a hole, Vp is given by: 




-^ V = -^ = [I pE 


(3.11) 




ITir 


where t is the mean free path between collisions and n is the 
electron or hole mobility. [Ref. 21, p.l7] 
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Whenever there is a carrier gradient or an electric 
field, current will exist in a semiconductor. The 
contribution of total particle current in the semiconductor is 
from both the holes and electrons flowing in opposite 

directions. The current density J due to an electric field 
(drift current of electrons) is then given by: 

J = e(n\i^ + piip)E (3.12) 

Since from Maxwell's equations: 

J = oE (3.13) 

where a is the conductivity, and can be expressed as a 
function of the material's characteristics: [Ref. 21, pp. 21- 
23] 

o = e(noii^ + Poji^) (3.14) 

When excess carriers are introduced in a semiconductor, 
the carrier density is no longer the equilibrium value. These 
excess carriers can be generated from external photoexcitation 
or any ionizing pulses incident on the material. This creates 
a carrier gradient or a local concentration of charge. These 
localized concentrations of charge follow Pick's law of 
diffusion which result in electron and hole particle diffusion 
current, and respectively: 

Jn = -e(-D„Vn) , fp = e{-DpVp) (3.15) 

where D^ p is the diffusion coefficient for electrons and holes 

in silicon. Then the total particle current due to diffusion 
and any electric field present is: 

Jp = e{n\ipE + DpVn + p[ipE - Dp Vp) (3.16) 

Utilizing charged particle transport theory, the kinetic 
theory of gases and the definition of mobility (Equation 
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3.11), the following Einstein relations can be derived: [Ref, 
21, pp. 22-23] 

3l = ^ = kT/e (3.17) 

\^n \^p 

Thermal carriers are produced when an electron in the 
valence band is thermally excited to the conduction band. 
This transition can take place in a single step (band-to-band) 
or through an intermediate band located in the forbidden band 
(this band is caused by the semiconductor impurities as 
discussed earlier). The latter process is most probable in 
silicon material. The rate of carrier generation can be seen 
to be proportional to the number of defects in the material. 
[Ref. 20, pp. 22-23] 

In order to maintain equilibrium carrier density in the 
material, a process called carrier recombination takes place. 
Like carrier generation, recombination is a two-step process 
in silicon like carrier combination. This process involves an 
electron in the conduction band combining with a hole in the 
valence band, hence destroying a hole-electron pair. A photon 
is released and in order to conserve both momentum and enercfy, 
a particle called a phonon is released (Refer to Figure 3.4). 
The phonon is analogous to the photon, but is associated with 
the acoustic waves of the vibrating lattice. Sometimes 
instead of photon emission, the energy is transferred to 
another electron as kinetic energy. This is known as an 
Auger process. In thermal equilibrium, the recombination rate 
and thermal generation rate are equal in order to maintain 
equilibrium charge density. [Ref. 21, pp.23-27] 

As mentioned before, the impurities in the semiconductor 
give rise to energy levels in the forbidden zone. This 
results in traps or recombination centers, where an electron 
can be captured or trapped and held by the impurity atom. If 
the energy level of the impurity atom has a high probability 


33 




(1) (2) 13) (41 

Electron Electron Hole Hole 

capture emission capture emission 


Figure 3.4. Carrier generation and 

recombination. From Ref[21, p. 24]. 


of reemitting the electron before it recombines with a hole, 
then this energy level is called a trap. This is the case if 
this impurity energy level is close to the conduction band. 
On the other hand, if the impurity energy level is midway in 
the conduction-valence energy gap, it is probable that the 
electron will recombine before it is reemitted. In this case 
it is termed a recombination center. [Ref. 21, p. 25] 

The average time an electron spends in the conduction 
band or hole in the valence band before recombination is 
called average or mean carrier lifetime, t . At equilibrium, 
it can be seen that: 




(3.18) 


where and are the average lifetimes for electrons and 
holes and and are the thermal equilibrium generation 
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rates. Under non-equilibrium states, generation and 
combination rates are not necessarily equal: [Ref. 21, p. 27] 


9p' 



(3.19) 


The concept of generation and recombination leads to 
definition of the carrier lifetime: "the average length of 
time that the charge carriers exist between the time that they 
are generated and the time that they recombine." [Ref. 23, p. 
105] The minority carrier lifetime, Tg and for electrons 

and holes respectively, can be expressed quantitatively by the 
following equation: 


T 


e 


1 

c N ' 

^ir^t 


1 


(3.20) 


where and Cp are the electron and hole capture 

crossections. [Ref. 19, p. 36-17] 

When the thermal equilibrium carrier densities are 
altered by some outside energy source, excess carriers are 
introduced into the semiconductor. This is called carrier 
injection. When the excess carrier concentration is less than 
that in thermal equilibrium, it is referred to as low-level 
excess carrier injection. High level carrier injection refers 
to carrier concentrations higher than the equilibrium 
concentration. [Ref. 22, p. 341] 

C. FORMATION AND CHARACTERISTICS OF A PN JUNCTION 


When n-type and p-type materials are joined together, 
electrons from the n-side diffuse to the p-side and holes from 
the p-side diffuse to the n-side. This forms what is called 
as a depletion or space charge region in the vicinity of the 
junction. There are positive ions on the n-side and negative 
ions on the p-side. This results in a built-in electric field 
across the region. The width of this depletion layer for an 
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applied voltage, is given by: 

1 = [2ee„((|)^ - V^) (3.21) 

where (|)o is the electrostatic potential between the two 

regions. Figure (3.5) illustrates the formation of the space 
depletion region. [Ref. 24, pp. 39-48] 
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Figure 3.5. Formation of the space charge 

(depletion) region in a pn junction. From Ref [21, 
p. 71] . 


Since there is a volume surrounded with opposite charge, 
a junction capacitance is generated. Because the capacitance, 
depends on the junction width, it is dependent on the applied 
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voltage, as in the following relationship: [Ref. 21, p.79] 


o \ 1/2 

^ 2 (1/JV^ + l/N^) ((|)^ - Vo) 


[3.22: 


When a junction is forward biased, the minority carriers 
are reorganized outside the depletion layer. This results in 
diffusion capacitance: 

Cd = + i';:-npo)exp(-^) (3.23) 

Diffusion capacitance can be related to minority carrier 
lifetime by: [Ref. 21, p. 80] 

= {e/kT)Ix (3.24) 

The bipolar junction transistor (BJT) is formed with two 
pn junctions as shown in Figure (3.6). The BJT consists of 
three parts: the base, emitter and collector. Generally the 
emitter is more heavily doped than the base to increase the 
transistor current gain, as will be shown later. For 
operation in the active region (for use as an amplifier) , the 
base-emitter junction is forward biased and the collector-base 
junction is reversed biased. 

In the active region of an npn BJT, the emitter-base 
current is primarily a result of the emitter injecting (or 
emitting) electrons into the base. The electrons, now 
minority carrier, diffuse through the base and flow into the 
reversed bias collector-base depletion region. Due to finite 
basewidth some electrons in the base will recombine with the 
holes in the base (assuming the basewidth is too narrow for 
recombination to occur) . The electrons flow into the 
collector (positive with respect to the base), where they 
appear as drift current. The collector current can be shown 
to be proportional to exp (eVs^/icT) and is given by: 
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= i'sexp ( eV^^/kT) 


(3.25) 


where Vg^ is the base-emitter voltage and is the saturation 

current. This shows that small fluctuations in the base 
voltage result in large fluctuations in the collector current, 
hence current amplification. This current amplification or 
gain is given as the ratio of the collector and base currents, 
P=I’e/-^s- Figure (3.6) shows the flow of currents in a npn BJT 
biased in the active mode. [Ref. 25, p. 196] 


Electrons 



Figure 3.6. Flow of currents in a NPN transistor 
biased for active operation. From Ref [18, p. 
339] . 


D. RADIATION EFFECTS IN A BIPOLAR TRANSISTOR 

The two basic effects of radiation on semiconductor 
material are injection of carriers due to ionization and 
displacement of atoms in the lattice structure [Ref. 2, p. 
1512]. The first of these, ionization, is most important for 
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this research, since with electron bombardment, ionization is 
the primary source of radiation damage. 

Ionization is the result of high energy photons 
interacting with the semiconductor and producing charged 
particles. The primary photocurrent will increase the base 
current resulting in a secondary photocurrent in the 
collector. The photoelectric effect, Compton scattering and 
pair production, described in Chapter II are included in the 
ionization processes discussed here. In turn, these charged 
particles generate photocurrents. Enhancement of bulk 
conductivity is a major effect of ionization and it is given 
by: [Ref. 21, p. 22] 


6a = e(|ij,6n + |ip6p) (3.26) 

where 6a is the difference between the post and pre-radiation 
values of conductivity, and bp,n are the excess carrier 
concentrations. 

Because BJT current gain is so important in analog 
networks (especially amplifier circuits), its change due to 
the radiation effect on it is of great concern. Current gain 
degradation can be the dominant effect in many bipolar 
transistors [Ref. 26, p. 2027] . A major cause of current gain 
degradation is the increase of surface states at the oxide- 
silicon interfaces. This in turn increases the surface 
recombination velocity, S^. Sp is a measure of the tendency 
of carriers to drift toward the surface and recombine. This 
increase in Sp results in a decrease in current gain P which 

can be seen from the following equation for p as a function 
of the transistor's physical characteristics: [Ref. 27, p. 
567] 
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(3.27) 


1/p = + o^W/o^L^ + A (P;r/i,^)2 

where: is the recombination velocity 

is the area for recombination 
JV is the base width 

is the base diffusion constant 
is the emitter area 
L is the diffusion length 

g is the conductivity of base and emitter 

The first term is due to surface recombination, the second is 
due to emitter efficiency and the third is due to volume 
recombination. 

Ionization can also cause high-injection level minority 
current into the base side of the base-collector region 
resulting in less exposed space charge from the base and more 
space charge on the collector side. This results in an 
expansion of the effective base width, and hence a decrease in 
P as reflected in the above equation. An increased base 
width also results in increased base transit time which 
decreases the transition frequency, fj.. [Ref. 21, p. 108] 
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IV. EXPERIMENTAL APPROACH 


A. OVERVIEW OF THE NPS LINAC 

This research utilized the 110 MeV NPS linear accelerator 
(LINAC) to simulate the total dose ionization effects on a 
device. The LINAC consists of three ten-foot accelerating 
sections, each powered by a 22 megawatt (peak) kylstron 
amplifier. An electron gun injects 80 keV electrons into the 
accelerator. The three klystrons modulate the electrons in 
"bunches" or packets. The klystrons are controlled by a 
master RF oscillator which turns on the beam voltage to the 
klystrons. The klystrons are driven by a tunable magnetron at 
2.856 GHz and at 6 kilowatts (peak). Since, at 80 keV the 
electrons are travelling about half the speed of light, those 
electrons not travelling at the same speed are lost. The 
number of electrons lost is reduced by inserting a 
"prebuncher" in the accelerator with a net result of reducing 
the special distribution of the electrons into a smaller 
bunch. The frequency of the electron beam is 60Hz with a 
pulse width of about 3.5 microseconds. Figure (4.1) gives an 
overview of the LINAC's major components and spatial layout. 

At the end of the third section of the accelerator is a 
collimator which provides a "source point" for the deflection 
magnet and quadrapole lenses. The deflection magnet separates 
the various energies in the beam, so that the desired energy 
can be selected by the energy defining slits. This magnet 
also serves to deflect, as the name implies, the electron beam 
so that the counting equipment is not in line with the 
undesirable forward radiation of the electrons. 

In order to measure the fluence (total number of 
electrons per unit area) of the electron beam, a secondary 
emission monitor (SEM) is utilized. This is located at the 
end of the target chamber of the LINAC. The electron beam 
passing through the SEM causes charge to accumulate on a 
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Figure 4.1. Spatial layout and major components of the 
NFS linear accelerator. From Ref. [28, p. 73]. 

capacitor which is remotely connected to a voltage integrator 
in the control room. The accumulated charge Q, due to N 
electrons of charge e, on the SEM is given by: 


Q = Ne = CV 


(4.1) 


[Ref. 29, pp. 20-23] [Ref. 30, pp. 7-8] [Ref. 31, p. 41-43] 

B. DOSIMETRY 

1. Overview 

Now that we have a radiation source, the problem arises 
on how to measure the radiation absorbed by the target device. 
There are many different methods of determining dose along 
with a variety of materials to accomplish this measurement. 
Among some types of dosimetry systems are: calorimeters. 







films, cobalt glass, Victoreen thimble-ionization chambers, 
thermoluminescent dosimeters (TLD's), to name a few. 
Selection of the right dosimetry system is dependent on the 
type of radiation testing done (i.e. ionizing, displacement 
damage), the radiation source and the target material (silicon 
in the case of this research). [Ref* 7, p. 87] 

Before going on to explain the method of dosimetry used 
in this research, a few terms relating to dose must be 
defined. First, is absorbed dose or total dose. This is the 
mean energy absorbed per unit mass of irradiated material and 
is given by D in the following equation: [Ref. 32, p. 1-3] 

D=LEo/^m (4.2) 

where AFTj, is the energy absorbed by the material. Dose is 

given in rads (radiation absorbed dose) , where one rad is 
defined as 100 ergs of energy absorbed per gram of material. 

Particle fluence, (j) is the number of particles, N 
incident on a cross-sectional area. Energy fluence, i|; is the 
radiant energy incident on a cross-sectional area. These two 
parameters are given by: [Ref. 32, p. 1-7] 

(|) = dN/dA, t|; = dE/dA (4.3) 

Another term of great importance in dosimetry, is charged 
particle equilibrium (CPE) . CPE is said to exist if the total 
energy carried out of a material by charged particles is equal 
to the energy carried into the material by charged particles. 
In the case of a photon beam incident on a material, CPE 
exists when the thickness of the material is approximately 
equal to the range of the most energetic secondary electron. 
This is sometimes referred to as the equilibrium thickness. 
Using this concept of CPE, the following equation can be 
derived: 


43 



( 4 . 4 ) 


D, ■ 


where (p-en/P^E.z the mass energy absorption coefficient at 
a photon fluence of energy E and material Z. is the 

absorbed dose of material Z. This equation assumes a 
monochromatic photon beam and little interference with 
surrounding material. [Ref. 33, p. 1171] 

Unfortunately, the concept of CPE is not valid in an 
electron beam environment. But a similar relationship can be 
derived for electron beams at high energies (E > 12 MeV) , 
assuming that electrons interacting with the'material will 
deposit energy within the material: 




(4.5) 


where ^col /P the mass collision stopping power expressed in 

MeV-cm^/g. Additionally, the assumption is made that the 
material is thin enough to ignore any bremsstrahlung 
absorption. [Ref. 33, p. 1172] 

2. Application of TLD Dosimetry 

In this research, manganese-doped calcium fluoride 
(CaF 2 :Mn) was used for dosimetry measurements. This is a 
synthetic material. Naturally occurring deposits of fluorite 
cannot be used because it displays non-reproducible light 
output and glow curves. This material is available in many 
physical forms and packaging. The type used here were 1/8 x 
1/8 X 0.035 inch chips. [Ref. 34, p. 3-9] 

When CaFj is exposed to radiation, holes and electrons 
are introduced and trapped by the metastable energy levels or 
traps in the forbidden gap (formed by the introduction of 
manganese). The density of the filled traps is proportional 
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to the total dose absorbed. To measure this absorbed dose, 
the TLD is then removed from the radiation source and heated 
in the TLD reader chamber. This heating gives the electrons 
in the traps enough energy to be released and combine with 
free holes, thus emitting light. The intensity of this light 
is proportional to the density of filled traps which in turn, 
is proportional to the total dose. [Ref. 32, p. 1-13] 

There are two major reasons for choosing this type of 
dosimetry. First, because of its linearity up to 100 krads, 
dosimetry calculations become very simple. Second, the energy 
dependence of unshielded CaFj is low at energies greater than 
one MeV. At energies less than this (not applicable in this 
research) , the TLD can be surrounded by a thick and thin 
metallic shield to flatten the energy dependence at energies 
below one MeV where error can be limited to 20% at energies 
greater than 80 keV. Figure (4.2) illustrates this energy 
dependence. [Ref. 34, p. 3-14] 

At high energies, the absorbed dose for CaF 2 can be 
considered equivalent to that of silicon. However, for 
energies less than about 12 MeV, Equation (3.4) must be used 
to find the equivalent silicon dose, given the CaF 2 dose. The 
energy absorption of TLD's and aluminum is plotted in Figure 
(4.3). Note that at 10 MeV (right boundary of the figure), 
the energy absorption for silicon and CaF 2 are almost 
identical. [Ref. 34, p. 1172] 

In conducting our dosimetry calibrations, the TLD chip 
was placed next to a "dummy" device under test (an identical 
device to the one actually irradiated) . This ensured that the 
calibration was done at the precise location of the device 
under test. The ideal situation would be to have the TLD 
placed inside the packaging material of the "dummy" chip and 
then pull it out after irradiation to measure the absorbed 
dose. This was obviously impractical for our setup, therefore 
the former method had to be used. 
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Energy, KeV 

Figure 4.2, Energy dependence of CaFjrMn. The top curve 
is the response for unshielded CaF 2 :Mn. The bottom curve 
is for shielded CaFjiMn. From Ref [34, p. 3-14] . 

Once the TLD was in place, the LINAC was pulsed (the beam 
was tuned and focused prior to this) for a predetermined 
amount of accumulated charge on the SEM, corresponding to a 
predetermined voltage on the capacitor. The TLD is removed 
and the total dose accumulated is measured. Several runs are 
made at various capacitor voltage levels. Since the SEM is 
known to be linear, the results of total dose vs. capacitor 
voltage can then plotted and extrapolated to get a value for 
Mrads per total voltage on the capacitor. Note that since the 
TLD's utilized, saturate at approximately 10 krads, the 



Energy Absorption In SlUoon Reitlive to the 
Energy Absorption In the Indlcetsd Materia} 



Figure 4.3. Energy dependence of various materials 
relative to silicon. From Ref [35, p. 1755] . 

irradiation time for the TLD's is on the order of seconds. 
According to simulations using CEPXS/ONELD radiation transport 
code to model the TLD response, the worst case deviation in 
total dose (Si) measurement may only be 15% at 60 MeV. This 
deviation is due to the differences in surrounding material in 
the case of the TLD chip and the silicon die. The actual 
experimental dosimetry plots are shown in the next chapter. 
[Ref. 33, p. 1174] 

C. DESCRIPTION OF THE DEVICE UNDER TEST 

The devices used for studying the total dose effects were 
Harris UHF-1 bipolar, 40 pin side brazed DIPs with gold plated 
kovar lids braze sealed. Each chip contained 13 NPN 
transistors of various emitter sizes. 
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The UHF-l process is a dielectric isolated, wafer bonding 
process that utilizes Silicon-On-Insulator (SOI) technology 
and vertical trenches. This technology reduces collector-base 
capacitance by forming a buried layer and contact and a 
recessed oxide region. The bonded wafer oxide allows for 
vertical inter-device isolation. Figure (4.4) reveals the 
cross section of the UHF-l architecture. Typical parameter 
values are shown for both the NPN and PNP transistors in Table 
4.1. [Ref. 36, p. 260-261] [Ref. 37, p. 37-42] 



Isolation Oxide 


Handle Wafer 

Figure 4.4. Cross section of the UHF-l architecture. 
From Ref [36, p. 261]. 

D. CIRCUIT DESIGN 

The 13 transistors were utilized for 3 separate sets of 
measurements. Nine of the transistors were utilized to 
construct a simple operation amplifier in a closed loop 
configuration. The vs. family of curves for different 
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PARAMETER 

DEFINITION 

NPN 

PNP 

UNITS 

Beta 

DC Current Gain 

100 

40 


Va 

Early Voltage 

60 

20 

V 

BVpBrt 


5.5 MIN 

4.5 BflN 

V 

BV,.„ 


12 MIN 

12 MIN 




8MIN 

12 MIN 

V 

Rirr 

Extrinsic Base Resistance 

40 

-25 

Ohms 

RpT 

Zero Bias Intrinsic Base Resistance 

350 

' 180 

Ohms 

R^ 

Parasitic Collector Resistance 

35 

50 

Ohms 


Intrinsic Collector Resistance 

225 

225 

Ohms 

Ctr 

Zero Bias Collector-Base Capacitance 

60 

80 

fF 


Collector-Substrate Capacitance 

30 

30 

■ fF 


Transition Frequency 

9 

5.5 

GHz 


Table 4.1. Typical parameters for UHF-1 transistor. From 
Ref. [36, p.263]. 

values of were monitored to extract current gain 

measurements of three transistors. The remaining transistor 
was used in a common emitter configuration for measuring the 
base and collector small signal currents as a function of 
frequency. These currents were used to calculate the 
transition frequency, fj.. 

1. Design of the Operational Amplifier 

The obstacle in designing an operational amplifier was 
the restriction to 9 NPN transistors and no PNP's. 
Additionally, the biasing of the transistors had to be 
calculated very carefully since these transistors had very low 
current ratings (the maximum current rating or power 
dissipation is not known). After some trial and error and 
destruction of transistors a suitable opamp design was found. 
Circuit simulation using PSpice was used extensively to model 
these transistors based on the specification sheet from Harris 
and experimental measurements. These transistor models in 
turn, were then used to design the opamp. The final design is 
shown in Figure (4.5) and its resulting phase plot and 
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frequency response in a closed loop configuration are shown in 
Figures (4.6) and (4.7). It was designed to have 
approximately a 1.5MHz gain-bandwidth product. 


R2> le.Sk 


R1 |20.7I< 




CIO t 


q3417 i <53417 


Figure 4.5. Operation amplifier design used for radiation 
testing. 


2. Common Emitter Circuit 

One way to experimentally determine the transistor 
transition frequency is to measure the ratio of the collector 
and the base small signal currents as a function of frequency. 
This was accomplished by building a common emitter circuit as 
shown in Figure (4.8). 

If Rg is large and is small compared to R^ then the 
collector and base currents are given as follows: 






Q P(V(VOUT^/ V(V4;+)) 

Frequency 


Figure 4.6. Phase plot of opamp. 


= 





V 



(4.6) 


A small signal of 10 mVp can be injected into this circuit at 
varying frequencies and the small signal base and currents can 
be measured. Then the AC P or P can be calculated from: 


= 



(4.7) 


and a plot for as a function of frequency can determine 
f^. The value of fj. can now be used to find the sum of the 
parasitic capacitances and in the high frequency hybrid- 
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3. Measurement of DC Current Gain ( P) 

A Tektronics 576 curve tracer was utilized to display the 
Ic vs. family of curves. The DC current gain, P was found 

for a given at a given ^CE by the following definition: 
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^CC 



Figure 4.8. Common emitter circuit used in experimental 
setup. From Ref. [25, p. 522] . 

6 = (4.9) 


While this measurement is very straightforward, the problem of 
noise had to be dealt with since the device and the curve 
tracer were connected together by over 80 ft of cable. The 
method used to deal with the noise is explained in the next 
section. 

E. EXPERIMENTAL SETUP 

The three circuits just described were implemented on a 
printed circuit board. With the exception of the common 
emitter circuit, all measurements had to be made in the LINAC 
control room which was approximately 80 ft away from the 
target device. This introduced an obvious noise problem. 
Since noise was a significant factor, all inputs and outputs 
to the circuit board were required to be connected with 
coaxial cables right up to the board connection. All the 
coaxial cables between the control room and the device were 
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triple shielded and ran through a steel conduit for most of 
the length. A setup of the circuit board is shown in Figure 
(4.9) . 

An AC small signal voltage (5-10 mVpp) was required for 
both the common emitter and the operational amplifier 
circuits. The signal source for both of these circuits was 
the output of the tracking generator of a Hewlett Packard 
3585B spectrum analyzer. This tracking generator had an 
option to manually sweep the frequency, which was ideal for 
the common emitter circuit measurement of and as a 

function of frequency. In the continuous sweep mode, the 
signal generator was used to display the closed loop frequency 
response of the opamp. 

The tracking generator output level was set to its 
maximum level of l.lVpp, to minimize noise interference with 
the signal. Of course it was necessary to step down the 
voltage for the opamp and common emitter circuits. This was 
done on the board utilizing the opamps to reduce the voltage 
to the 10 mVpp range. This worked extremely well as can be 
seen by the closed loop response of the opamp in the next 
chapter. This is the display that was present while the LINAC 
klystrons were on (the klystrons were the most significant 
noise source in the LINAC environment). An oscilloscope was 
used to remotely measure the opamp slew rate. 

Because of the small amplitude of the voltages that were 
measured in the common emitter circuit, auto-ranging 
voltmeters were setup near the target device. The voltmeters 
were read remotely in the control room via closed circuit 
video cameras. Figure (4.10) illustrates this setup. 

As stated previously, two Tektronics curve tracers were 
utilized to obtain current gain (refer to Figure (4.11)). 
With some experimentation it was found that negligible noise 
was observed on the cu 2 rve tracer display if the inner 
conductor of the coaxial cable was used for the stepped base 
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current. The emitter (grounded at the curve tracer) was 
placed on the outer conductor of the coaxial cable. The 
collector current was place on a single conductor cable. This 
arrangement produced the optimum display on the CRT. There 
was negligible difference between the remote curve tracer 
measurements and local curve tracer measurements. The quality 
of the traces can be seen in Figure (4.12). 
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Figure 4.9. Setup of circuit board in the target area of 
the LINAC. 












































Figure 4.11. Two Tektronic 576 curve tracers utilized for 
measuring current gain during irradiation. 
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Figure 4.12. Displays of Tektronic 576 curve tracer while 
irradiating device-under-test. 
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V. RESULTS 


When analyzing the experimental results, it must be kept 
in mind that only the transistors in the circuit were 
irradiated and the other components were off chip and 
discrete. In the case of the opamp, the compensating 
capacitor and resistors received an unmeasured, but 
insignificant dose. This research concentrated on degradation 
of opamps due only to radiation damage in the transistors. 

A total of five days of irradiation runs were conducted, 
with two days using 30 MeV electrons and three days using 60 
MeV electrons. The reason for varying the electron energy was 
twofold. First, we wanted to see if silicon devices had any 
energy dependence in the 3 0 MeV and 6 0 MeV range. Secondly, 
at 30 MeV, it was not possible to achieve a high enough dose 
rate (without bringing radiation-free areas above 
regulations), in order to degrade the transistor P 
appreciably in the time allotted. At 60 MeV, the beam is more 
focused, hence less electrons leak out and contaminate the 
radiation-free areas. 

The first run on 2 February, chip #34 was irradiated with 
30 MeV electrons in the circuit configurations specified in 
Chapter IV (as all the runs were) . The results of the 
dosimetry are shown in Figure (5.1). The total dose 
accumulated in the silicon chip was 4.7 Mrads(Si) . The effect 
of radiation on the current gain, P (h^) is shown as a 
function of total dose for three of the 13 transistors in 
Figure (5.2). Note that P is normalized to its initial 
value, Po- Total dose in rad(CaF 2 ) is taken to be equivalent 

to total dose rad (Si) (Refer to Chapter IV) . The current gain 
decreased by an average of 15% for the three transistors. 
From previous research that modeled radiation effects on 
opamps, we didn't expect any measurable degradation in the 
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Figure 5.2. Closed loop response of opaitip 
after transistors exposed to 4.7 Mrads(Si). 
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Figure 5.3. Current gain degradation as a 
function of total dose. HFE^ is the value at 
the start of the trial. 

opamp parameters. Figure (5.3) shows the closed loop response 
for the opamp with an accumulated dose of 4 Mrad(Si). This 
response was monitored during the entire time of irradiation 
with no noticeable changes in either the 3dB frequency, 
Bandwidth (BW) or gain. The slew rate was checked 
periodically, and as expected, there was no change. 

It is interesting to compare the effects of different sources 
of ionizing radiation. Figure (5.4) compares our experimental 
data from Figure (5.3) to experimental data of same type of 
chips irradiated with a Shepherd Co60 source. The data points 
label #13HF47_83 and #UHF47_84 were the transistors irradiated 
with the Co60 source up to a total dose of 1 Mrad(Si) [Ref. 
38]. Of course, all we can say from this comparison is that 
both radiation sources have similar effects on these type of 
devices. Obviously, more extensive testing with controlled 
conditions would be required to make any kind of quantitative 
relationship from experimental data. 
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Figure 5.4. Comparison of effects of 
irradiation of UHF-1 transistors by a Co60 
source and LINAC. 


For the next radiation testing, it was decided to 
continue the testing again at 30 MeV with the previously 
irradiated chip #34, to see if we could bring the current gain 
down to 20% of its original value. This time the total dose 
achieved for the day was 17.1 Mrad(Si) , a two day combined 
total of 21.8 Mrad{Si). The dosimetiry results for the 
remaining trials are shown in Figure (5.5). 

As in the previous trial, there was no observable change 
in the parameters of the opamp. Figure (5.6) represents the 
current degradation relative to the start of this particular 
run. The total dose scale represents the cumulative total 
dose, including that of the previous day. Figure (5.7) 
combines the data from the two days of radiation exposure as 
if it were taken in one continuous trial. At the interface of 
the two trials, a discontinuity can be seen in the data. It 
would appear if some annealing took place overnight 
(approximately 14 hours between trials). However, this is 
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6 Feb Dosimetry 


7 Feb Dosimetry 



mV 


mV 


Figure 5.5. Dosimetry results for 6 through 9 February 
trials. 

probably due more to thermal action with a little contribution 
to annealing. Temperature was not monitored in this 
experiment, but probably should be in any future testing. 
Previous testing with these devices suggest that little 







Figure 5.6. Current gain degradation as a 
function of total dose. HFE,, is the value 
at the start of the second trial. 



Total Dose (Krad (CaF2) ) ^ 

Figure 5.7. Results of two different days 
of irradiation with 30MeV electrons. HFE^ is 
the value at start of the first day. 
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annealing does take place. The transistors were still between 
50 and 70 percent of their original gain after a two day total 
dose of 21.8 Mrad(Si) . This still was not enough degradation 
to see a change in the operational amplifier characteristics. 
The decision was made to irradiate chip #34 a third day to 
degrade the current gain even more. In order to achieve a 
higher dose rate, ultimately a higher total dose in the time 
allotted, it was decided to increase the electron beam to 60 
MeV. As alluded to earlier, increasing the particle fluence 
(to achieve a higher dose rate) at 30 MeV increased radiation 
levels to unacceptable values outside the laboratory. 

Due to a damage to a transistor in the common emitter 
(CE) circuit, transistor #34_10 was used as a replacement in 
the CE circuit. Therefore the remaining plots will have only 
two transistors plotted for 

The results of the 7 February trial are plotted in Figure 
(6.8). The three day total dose accumulation was 163 
Mrad(Si). This third trial alone contributed 141.2 Mrad(Si), 
utilizing 60 MeV electrons. As before. Figure (5.8) reflects 
a current gain normalized to the starting value on the third 
day, while the total dose scale reflects the three day total. 
Figure (5.9) plots data for the trials as if it were one 
continuous trial. It can be seen that there was no 
"annealing" effect at the boundary between the second and the 
third trials. From this plot, it appears that the 60 MeV 
electrons did not change the trend of the current degradation. 
So far the current gain of the two transistors are still about 
30% of their original value after receiving 160 Mrad(Si). 

There was still no effect observed on the corresponding 
opamp characteristics. The gain, slewrate and 3dB frequency 
have not changed. Testing proceeded on new chip #30 in the 
same manner as before, but this time only 60 MeV electrons 
were to be used. 
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Figure 5.8. Current gain degradation as a 
function of total dose. HFEg is the value 
at the start of the trial. 



Figure 5.9. Results of three different 
trials utilizing chip #34. 
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The results of the first run on 8 February with chip #30 
are plotted in Figure (5.10) . The total absorbed dose at the 
end of this trial was 202 Mrad{Si). The current gain of the 
transistors was decreased to 25% and 35% of the original value 
for #30_6 and #30_9 respectively. As before, this still 
wasn't enough to observe any degradation in the opamp 
parameters. Therefore a second trial at 60 MeV was conducted 
on 9 February to achieve a combined two day total dose of 500 
Mrads{Si). As before, these results are plotted in Figure 
(5.11) just for the second trial, and for the combined two 
trials, in Figure (5.12). 

It can be seen that up until approximately 60 Mrad(Si), 
the two transistors degrade at almost exactly an identical 
rate. After about 60 iy[rad(Si) total dose, the rates diverge 
significantly. Referring back to Figure (5.9), it can be seen 
that Chip #34 behaved almost the same way, at the same total 
dose. Chip #30_9 has an emitter perimeter-to-area (P/A) ratio 
that is 6 times that of chip #30_6, therefore one would expect 
that chip #30_6 and chip #34_6 to degrade faster than 
#34,3 0_9, according to Equation (3.27) . Note that chip #30_6 
and #34_6 have the same P/A ratio as do chip #30_9 and #34_9. 

Another observation can be made by comparing the results 
from two different energies of incident electron beams. Since 
the first run irradiated chip #34 to 4.7Mrads(Si) continuously 
at 30 MeV, this set of data will be plotted along with that of 
chip #30 which was irradiated continuously at 60 MeV to 202 
Mrads(Si). The two transistors of each chip correspond to 
each other, therefore the P/A ratios are the same for each 
counterpart. This is shown in Figure (5.13). 

In order to compare the effects of electron beam energy 
and the effects of removing the device from the radiation 
source for for over 16 hours. Figure (5.14) is plotted. In 
Figure (5.14) , the current gain data is plotted for both chips 
up to 160 Mrads(Si). It is important to note that the plot 
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Figure 5.10. Current gain degradation as a 
function of total dose. HFE„ is the value of 
the start of the trial. 



0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 


dose (Krad(CaF2)) ^ 

Figure 5.11. Results two different trials 
utilizing chip #30. 
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Figure 5.12. Current gain of both irradiated 
chips as a function of total dose up to 4.3 
Mrads (Si). 



Figure 5.13. Current gain degradation as a 
function of total dose. HFEq is the value at 
the start of the second trial. 




















for both transistors on chip #30 were irradiated continuously 
in a 60 MeV electron beam for the entire range of total dose 
plotted. On the other hand, the transistors on chip #34 were 
irradiated in three separate steps: using a 30 MeV electron 
beam up to 4.7 Mrads(Si) , using a 30 MeV beam with the initial 
4.7 Mrads(Si) and continuing to a cumulative total of 21.7 
Mrads(Si) , and finally using a 60 MeV beam to bring the final 
cumulative dose to 160 Mrads. The intervals between each step 
offered the device a "rest" time of 18 to 24 hours without 
radiation exposure. 

Nothing conclusive can be drawn from this plot, however, 
one can see that the 3-step plot of the chip #34 transistors 
almost appears to be continuous. It does show a sharper rate 
of gain degradation than chip #34 transistors, but this could 
be due to the fact that they are two distinct chips which are 
bound to have differences even if they were made on the same 
wafer. This can be verified by further irradiation of the 
same types of devices without any "rest" times from radiation 
exposure. 

Although not conclusive, there is strong evidence from 
these trials that current gain degradation vs. total dose 
relationship is not energy dependent in the 30 to 60 MeV 
range. Again, strong evidence to support this statement can 
be accomplished by irradiating continuously the same type of 
devices under a 30 MeV beam and another at 60 MeV. 

Not much can be said about the opamp other than no 
measurable change was observed in either chip #3 0 or chip 
#34's opamp circuit. This is not completely unexpected, for 
it supports previous thesis research conducted on modeling 
radiation effects with PSpice. Figure (5.15) is the result of 
a simulation that degrades the current gain (or P in the 
figure) down to several levels or percentages of remaining P . 
A simulation is run for a generic ^A741 opamp utilizing these 
degraded transistor models. It can be seen that to obtain 
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dose (Krad<CaF2)) ^ 

Figure 5.14. Cumulative current gain data 
for both chips up to 160 Mrad{Si). 

Gain (dB) 

5 


Frequency (Hz) 


Figure 5.15. piA741 operational amplifier 
degradation curves as function of percent 
remaining current gain. From Ref. [3 9, p. 
65] . 
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observable degradation in the opamp gain, the transistor 
current gains must degrade to 5% of their original value. In 
this research, we were only able to bring the current gains 
down to an average of 30% of the original value, therefore we 
did not see any degradation in opamp parameters. 

It should be noted that we did not include data from the 
common emitter circuit, specifically the ac voltage 
measurements corresponding to the base and collector currents. 
This is because the data was not consistent. There were 
severe fluctuations when the klystrons were turned on and the 
initial calculations gave a value of that was many orders 
of magnitude lower than from the manufacturer's specification. 
There was indication that fj. did decrease with total dose as 

expected, but no useful data was obtainable. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


This research was intended to help create a more accurate 
model of bipolar transistors using data derived from 
experiment. Because of limited time and limited number of 
devices, the quantitative data obtained from this research is 
not useful for forming an accurate model. Many repeated runs 
under varying conditions and configurations would be required 
to form an accurate quantitative description. However, the 
procedures investigated can be used in establishing a baseline 
for further improved experiments. 

Some interesting conclusions can be drawn from this 
investigation. First, there is strong evidence that the 
current degradation in the UHF-1 Harris bipolar transistor is 
not dependent on the energy of the incident electron beam, at 
least between the energies of 30 MeV and 60 MeV. Previous 
research at NPS and the Jet Propulsion Laboratory has shown 
this to be true for other semiconductor devices in the energy 
ranges of 16-84 MeV [Ref. 40, p. 4] [Ref. 41, p. 11] [Ref. 42, 
p. 19]. Note that shielding effects are not included here. 

Secondly, this investigation showed that in-situ 
measurements for transistors and opamps are possible, 
utilizing a single die, thus ensuring a uniform radiation 
distribution. The resulting displays on the spectrum analyzer 
and curve tracers showed almost non-observable levels of 
noise, even with the LINAC kylstrons turned on. 

It was surprising that the transistors retained as much 
as 25-30% of their original gain after being exposed to over 
500 Mrads (Si). Similar work done using the NPS LINAC, found 
that integrated opamps of different technology survived after 
receiving doses of over 68 Mrads{Si) [Ref. 42, p. 92]. 
Survival under these high doses is questionable and should be 
further investigated with the same type of devices used in 
this research. Of interest is the data in Figure (5.4), where 
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the total dose response using a COgo source compared favorably 
with that of the NPS LINAC. However, nothing is known of the 
effects of 500 Mrad{Si) using a COgg source. 

Continuation of this research should include 
investigation into a noise resistant method to measure the 
transition frequency f^. This is necessary in order to 
understand the effect of radiation on the parasitic 
capacitances inherent in the transistor. In addition, further 
irradiation trials should include testing these SOI bipolar 
devices to the same doses used in this experiment, but without 
any "rest time" to see if this has an effect on the overall 
total dose response. 

In order to understand the effects of capacitances and 
resistors used in integrated circuits, these devices should be 
irradiated in a similar fashion as was done here. The ideal 
test device required to accomplish the goal of obtaining an 
effective model for an opamp would be a die with a standard 
opamp built on it along with individual transistors having the 
same dimensions and design as those used in the opamp. This 
will ensure that all components of the opamp and the 
individual transistors will be exposed to the same radiation 
and the responses to both can be compared. 
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APPENDIX. GLOSSARY OF SYMBOLS AND TERMS 


Due to differences in symbology and terms used in various 
sources, this appendix provides a compact list of symbols and 
terms along with their definitions as used in this paper. 

E - energy 

e - charge of one electron equivalent to 1.6X10'^® Columbs 

m - mass 

h - Plank's constant, 4.136x10'^® eV* sec 

V - frequency 

V - velocity 

Vq - stopping potential 

p - momentum 

A - wavelength, or magnetic (invariant) latitude, in the 

case of Mcllain's L parameter 
c - speed of light, 2.9979 x 10® meters/sec 

k - Boltzmann constant, 8.617 x 10'® eV/°K 

T - temperature (300K° usually considered room 

temperature) 

Eq - permittivity of free space, 8.854 x 10'^^ C^*N‘^'m’^ 

Z - atomic number 

n - neutron 

p - proton 

L - Mcllain's L parameter, approximately the giocentric 

distance of a field line in the geomagnetic equator 
R - radius of the earth 

M - geomagnetic dipole moment 

B - magnetic field 

Ef - Fermi energy 

Ng - density of states per unit volume 


77 



effective mass of an electron in a given state 


m 


* 

e 


irir 


no 

Po 

ni 
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'n,p 


9o 



1 
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Cd 


C.. 


a 


P 






Q 


- effective mass of a hole in in a given state 

number density of states for electrons (carriers/cmj) 

number density of states for holes (carriers/cm^) 

intrinsic carrier density 

donor concentration 

acceptor concentration 

mean free path between collisons 

carrier mobility 

current density, Amperes/m^ 

particle current 

electron, hole diffusion coefficient 
carrrier thermal equilibrium generation rate 
electron, hole capture cross-sections 
density of recombination centers (traps) 
width of depletion layer 
diffusion length 
junction capacitance 

diffusion capacitance 
collector-base capacitance 
emitter-base capacitance 
coductivity 

excess carrier concentrations 
BJT current gain 
recombination velocity 

base width 

transition frequency 
total electric charge 


78 



D 

t 

^col 

V-eJ9 


total or absorbed dose 
energy fluence 

mass collision stopping power 
mass energy absorption coefficient 
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